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We consider the evolution of a system composed of non-interacting point particles of mass 
m in a cylindrical container divided into two regions by a movable adiabatic wall (the adiabatic 
piston). We study the thermodynamic limit for the piston where the area A of the cross-section, 
the mass M of the piston, and the number of particles go to infinity keeping A/M and N/M 
fixed. The length of the container is a fixed parameter which can be either finite or infinite. In this 
thermodynamic limit we show that the motion of the piston is deterministic and the evolution is 
adiabatic. Moreover if the length of the container is infinite, we show that the piston evolves toward 
a stationary state with velocity approximately proportional to the pressure difference. If the length 
of the container is finite, introducing a simplifying assumption we show that the system evolves 
with either weak or strong damping toward a well-defined state of mechanical equilibrium where 
the pressures are the same, but the temperatures different. Numerical simulations are presented to 
illustrate possible evolutions and to check the validity of the assumption. 
Key Words: Liouville equation, adiabatic, piston, equilibrium, damping. 

1. INTRODUCTION 

The adiabatic piston problem is a well-known controversial example of thermodynamics where the two principles of 
thermostatics (conservation of energy and maximum of entropy) are not sufficient to obtain the final state to which an 
isolated system will evolve. This can be understood since the final state will depend upon the values of the "friction 
coefficients" , which however do not appear in the entropy function. Similarly, taking a microscopical approach, i.e. 
statistical mechanics, this implies that one can not find the (thermodynamical) equilibrium state by phase space 
arguments. 
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Let us recall the problem. An isolated, rigid cylinder is filled with two ideal gases separated by a movable adiabatic 
rigid wall (the piston). Initially the piston is fixed by a brake at some position Xq and the two gases are in thermal 
equilibrium characterized by their respective pressures and temperatures (or equivalently their energies ) . 
At a certain time, the brake is released and the question is to find the final equilibrium state assuming that no friction 
is involved in the microscopical dynamics (see [1-3] and references therein for different discussions). Experimentally 
the "adiabatic piston" has been used already before 1940 to measure the ratio of the specific heat of gases 7 — Cp/c^. 

Recently, the problem was investigated using the following very simple microscopic model [5-9] , which is a variation 
of the model already introduced in [p^ . The system consists of N non- interacting point particles of mass m in an 
adiabatic, rigid cylinder of length 2L and cross-section A. The cylinder is divided in two compartments, containing 
respectively and iV+ particles, by a movable wall of mass M 3> m, with no internal degrees of freedom (i.e. an 
adiabatic piston). This piston is constrained to move without friction along the x-axis. The particles make purely 
elastic collisions on the boundaries of the cylinder and on the piston, i.e. if v and V denote the x-component of the 
velocities of a particle and the piston before a collision, then under collision on the piston: 

v^v' = 2V-v + a{v-V) V = V + a{v-V) (1) 

where: 

2m 

a= (2) 

M + m ^ ' 

Similarly, for a collision of a particle with the boundaries at x = ±L, we have: 

V —> v' — —V (3) 

Since there is no coupling with the transverse degrees of freedom, one can assume that all probability distributions 
are independent of the tranverse coordinates. We are thus led to a formally one-dimensional problem (except for 
normalizations) and therefore in the model we shall assume that all the particles have velocities parallel to the x-axis. 

Since in physical situations m ^ M, this model was investigated in |^] in the limit where e = \frnjM <C 1 with M 
fixed. It was shown that if initially the two gases have the same pressure but different temperatures, i.e. p+ = pT 
but ^ , then for the infinite cylinder {L — 00) and in the limit e — )■ {M fixed, finite), the stationary solution 
of Boltzmann equation describes an equilibrium state and the velocity distribution for the piston is Maxwellian, with 



temperature Tp = T+, i.e. the piston is an adiabatic wall. However, to first order in e, the stationary solution 
of Boltzmann equation is no longer Maxwellian and describes a non-equilibrium state where the piston moves with 



constant average velocity V — ^jTimkB /?>{\/T+ — \/T^)/M towards the high temperature domain, although the 
pressures are equal. In it was shown that for L = oo, the evolution towards the stationary state is described to 
order zero in e by a Fokker-Planck equation, which can then be used to study the evolution to higher order in e. The 
case of a finite cylinder was investigated in |^] by qualitative arguments and numerical simulations. It was shown that 
the evolution take place in two stages with very different time scales. In the first stage, the evolution is adiabatic and 
proceeds rather rapidly, with or without oscillations, to a state of "mechanical equilibrium" where the pressures are 
equal but the temperatures different. In the second stage the evolution takes place on a time scale several orders of 
magnitude larger (if e <^ 1); in this second stage, the piston drifts very slowly towards the high temperature domain, 
the pressures of both gases remain approximately constant, and the temperatures vary very slowly to reach a final 
equilibrium state where densities, pressures and temperatures of the two gases are the same. It was thus concluded 
in that a wall which is adiabatic when fixed, becomes heat-conducting under the stochastic motion. However for 
real systems, the time involved to reach the thermal equilibrium will be several million times the age of universe and 
thus for all practical purposes the piston is in fact an adiabatic wall. 

On the other hand, for systems with e « 1, or e = 1 [II-I4], the system evolves directly toward thermal equilibrium 
and can never be considered as an adiabatic wall. 

Another limit has been recently considered. In p5[ |, the authors have studied the case where the mass m of the 
particles is fixed, but L tends to infinity together with M ^ L^, N ^ and the time is scaled with t ~ tL. They 
have shown that in this limit, the motion of the piston and the one-particle distribution of the gas satisfy autonomous 
coupled equations. 

Following a suggestion of J.L. Lebowitz jl^, we shall analyse in the following still another limiting procedure. We 
consider the thermodynamic limit for the piston where m and L are fixed, but the area A of the cross-section of the 
cylinder tends to infinity while 

2mA + 
7 — TT and R — m — — - (4) 

' M+m M ^ ' 

are kept constant ("— " refers to the left and "-f" to the right of the piston). 

The object of this article is to show that in the thermodynamic limit for the piston, the motion of the piston is 
"adiabatic" in the sense that it is deterministic, i.e. < V^" >t—< V >", no heat transfer is involved, the entropy of 
both gases increases, there is factorization of the joint distribution for one particle and the piston, and the system 
evolves toward a state of mechanical equilibrium, which is not a state of thermal equilibrium. Furthermore numerical 
simulations indicate that the evolution is such that the energy of the gas increases under compression (because work 



is done on the gas); moreover the simulations also indicate that the evolution depends strongly on for small values 
{R^ < 10) but tends to be independent of for larger values (see Figure 2 and 3). 

In section 2 we derive coupled equations for the one-particle velocity distributions. The thermodynamic limit for the 
piston is investigated in section 3. In section 4 we discuss the case where the length of the container is infinite; then 
the case of a finite container is considered in section 5. We present numerical simulations in section 6, and finally the 
conclusions in the last section. 

2. COUPLED EQUATIONS FOR THE ONE-PARTICLE VELOCITY DISTRIBUTIONS 

A. Initial condition 

We label (— ) those properties associated with the particles in the left compartment and (+) those associated with 
the right compartment. The number N~ and of particles in each compartment is fixed. For the sake of clarity, 
when necessary we denote {xi,Vi), i = 1 . . ■N~ , the position and velocity of the particles in the left compartment 
("left particles"), by {yj,Wj), j = 1 . . . iV+, the position and velocity of the right particles, and by (X, V) the position 
and velocity of the piston. Taking the boundaries of the cylinder at a; = ±L, we thus have: 

-L<Xi<X<yj<L (5) 

In this section, L, and M are finite. Initially, the piston is fixed at {X = Xa, V = 0) and the particles on both 
sides are in thermal equilibrium at respective temperatures T~ and T+, i.e. the probability distribution / in the 
whole phase space is given by : 

f{xi,vi; . . .■,XN-,VN-;X,V;yi,wi; . . .■,yN+,WN+;t = 0) = 

= [^(X + XoV" [^(£-Xo)]^+ n ri^^n^o - Xi)e{xi +L)fl f+{wj)e{y^ - X,)e{L - yj)S{X - X,)S{V) 

z=l J— 1 

(6) 

where 9 denotes the Heaviside step function. The probability distributions f^{v) are functions of v"^ with: 

/oo pO poo 

f'^{v)dv = l, 2m v^f+{v)dv = kBT+, 2m f {v)dv = UbT' (7) 

-cx) J — oo Jo 

For example, one could take for (u) Maxwellian distributions with temperatures T"^ , or functions which are non-zero 
only for \v\ G [vmin, Vmax]- For t>0, the piston moves freely and the problem is to study its evolution. 
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B. Liouville equation 



The distribution function / = f{{xi, Vi}; X, V] {yj, Wj}; t) at time t is a measure over the whole phase space, which 
is symmetric in the j-coordinates (resp. in the j-coordinates), with the foUowing normahzation which reflects the 
underlying dimension d = 3 of the system: 

(8) 

■' »=i i=i 

(with A the cross-section of the cylinder). Its evolution in time is given by the Liouville equation: 

1=1 j=i "■'J 

together with the boundary conditions: 

= =/(---)|..=-i,-., (10) 

/(■ ■ ■)\yj=L,Wj = /(• • ■)\yj=L~Wj (11) 

The right-hand side T of Eq. ^ takes into account the elastic collisions of the particles with the wall at ±L and with 

the piston. It decomposes into four terms. 

1 - Elastic collisions of left particles with the wall at x = —L: 

1i = Y.f=i ^{^i + L)vi[0(vi)f{xi,vi; ...■,Xi, -Vi] . . . ; x^- , Wjv- ; X,V;yi,wi; . . . ; y^+ , wn+ ; t) 

+ d{-Vi)fixi,vi; ...■,Xi,vf,...; x^- , v^- ; X,V;yi,wi; . . . ; yN+,WN+;t)] (12) 

Using the boundary condition (p^, this contribution is simply: 

AT" 



I^=J2S{x, + L)v,f (13) 

i=l 

2 - Elastic collisions of right particles with the wall at x ~ L: 
l2 = - J2f=i Siyj-L)wj[e{-'Wj)f{xi,vi]...;xN~,VN-;X,V;yi,wi;...;yj,-Wj;...; ?/jv+ , Wn+ ; t) 

+ 0{wj)f {xi,vi; . . . ; x^- , v^- ; X,V;yi,wi; . . . ■,yj,Wj; . . . ; yN+,WN+ ; i)] (14) 
and using the boundary condition (pi]), we have: 

l2 = -Y.6{y,^L)w,f (15) 



3 - Elastic collisions of left particles with the piston at x = X : 



^3 = Ljli ^i^i ~ ~ Vi)[0{V ~ Vi)f{xi,vi; . . ■■,Xi,v'^; . . . ; xn- ,vn-; X,V' ■,yi,wi; . . . ■,yN+ ,WN+-:t) 

+ e(vi - V)f{xi,vi; ...■,Xi,Vi;.. . ; xn- ,vn-; X,V;yi,wi; . . . ■,y^+ ,WN+;t)] (16) 
where (?;■, V) are given by Eqs. (0, 2). 

4- Elastic collisions of right particles with the piston at y = X : 

Ii=- ^iVj - X){V - Wj)[0{wj - V)f{xi,vi; . . .■,X]sr-,Vf^-;X,V';yi,Wi; . . . ■,yj ,Wj; . . . ;yf^+ ,WN+;t) 

+ 9{V - Wj)f {xi,vi; x^- , vm~ ; X,V;yi,wi; . . . ;yj,Wj; . . . ; yN+,WM+ ; t)] (17) 

The correlation functions for k left particles, / right particles (and possibly the piston) are defined by integration of 
the distribution / over (N^ — k) left particles, {N^ — I) right particles, and the piston coordinates (possibly not), 
multiplied by the factor: 

{N- -k)l {N+ -ly. ^^^^ 
The evolution for the correlation functions is then obtained from Liouville equation by integration over the corre- 
sponding coordinates. 

C. Equation for p^{x,v;t) 

We denote by p~(x,v;t) the single-particle distribution in the left compartment, obtained by integrating over all 
variables except xi = x and vi = v, with the normalization: 

N~ 

p {x,v;t)dxdv = (19) 

Its evolution involves pi.p the correlation fmiction for one left particle and the piston, with: 

pi^p{x,v;X,V;t)dXdV = p-{x,v;t) (20) 
Pi,p{x, v; X, V; t)dxdv = — ^'(X, V- 1) (21) 

where ^'(X, V] t) is the normalized probability distribution for the piston, obtained by integrating / over all variables 
except {X^V). Integrating Liouville equation (^ over all variables except xi — x and vi — v, together with the 
boundary conditions (|l^) yields: 



with the initial condition: 



{dt + vdx)p {x,v;t) = 5(x + L/2)vp {x,v;t) 

+ {V-v)[e{V-v)pi^p{x,v';x,V';t) + e{v-V)piA^,v;x,V;t)]dV (22) 



TV" 

p-{x,v;t^O)^f-{v)j^^^—^e{Xo~x)eix + L) (23) 



D. Equation for p'^{y,w,t) 

Similarly, we introduce p'^iy, w; t) the single-particle distribution in the right compartment and pp i the correlation 
function for one right-particle and the piston: 

N+ 

p+{y,w]t)dydw = — (24) 

pp,i (X, V;y,w; t)dXdV = p+ [y, w; t) (25) 

N+ 

pPsiX, V; y, w; t)dydw = — ^'(X, t) (26) 
Integrating Liouville equation (^ together with the boundary condition ( pl| ) yields: 
{dt + wdy)p+{y,w;t) = -S{y ~ L)wp+ {y,w;t) 

[V - wMw - V)pp,i{y, V- y, w';t) + 9iV - w)pp,iiy, V; y, w; t)]dV (27) 
together with the initial condition: 

N+ 

p+{y^w-t = Q) = f+{w)^^^——^e{y-X^)e{L-y) (28) 

E. Equation for l^;t) and ^{V;t) 

Integrating Liouville equation (^ over all variables except {X, V), and assuming ^![X — ±L, V; t) — 0, yields: 

/CX3 
{V - v)[e{V - v)pi,p{X, v'; X, V'; t) + 9{v - V)pup{X, v; X, V; t)]dv 
-oo 

/oo 
{V-w)[0{w - V)ppA{X,V';X,w';t) + e{V -w)pp,i{X,V;X,w;t)]dw (29) 
-oo 

together with the initial condition: 

'i'{X,V;t = 0) = S{X - Xo)d{V) (30) 



Finally integrating ( p9| ) over X leads to the evolution equation of the distribution function ^{V; t) for the velocity of 
the piston: 



dMV-t) =A {V- vMV - v)p-^^^{v'- V'-t) + e{v ~ V)p-^,f{v- V-t)]dv 

J —OQ 

-A I {V-wmw-V)pt^^f{w'-y-t) + e{V-w)p-t^^j{w-V-t)]dw (31) 

where: 

/oo 
p^AX.v;X,V-t)dX (32) 
-OO 

Pturfi^-^V-'t) = pp^,{X,V;X,w;t)dX (33) 



represent the joint distribution for the piston velocity V and the particle velocity v on the left (resp. right) surface of 
the piston. Let us note that 



/oo /-oo 
/ ptrA^;V;t)dvdV (34) 
-C30 J —CO 

represents the density of particles on the left (resp. right) surface of the piston. 

Let us transform equation (^ij) by considering its action < <J>, g > (t) = J^^^{V;t)g{V)dV on a test-function 
g{V). We have by definition: 

d 



^^<$,g>=y dMV;t)g{V)dV (35) 

Introducing the expression ( |3l] ) of dt^{V;t) in Eq. ( |35| ) and making the change of integration variables {v,V) 
{v',V') in the term involving /07tir/(^'' ^) ^^'^j respectively, the change of integration variables {V,w) {V',w') 
in the term involving p^urfi^''^ ^)); obtain: 

j^<<i>,g>=Aj dVdv(v ~ V)e{v ~ V)p-^^j{v- V- 1) [g{V + a{v-V))- g{V)] 

/OO 
dVdw{w - V)9{V - w)p+^^f{w; V;t)[g{V + a{w - V)) - g{V)] (36) 
-OC 

We then expand g{V + a{v — V)) in powers of a: 

g{V + a{v^V))- g{V) = ^li^—^^Cfc+i) (y) (37) 

and transform the terms as follows: 

aAj^^ dVdv{v - V)e{v - V)p;^^Jv; V;t){v - y)fe+ig(fe+i) (F) = 



/ 



("l^'^'nTnTfeTT / i^~Vr+'pJ^^^{v;V;t)dv ,g) (38) 
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with 



7 = aA ■■ 



2mA 
M + m 



In conclusion, the equation for the distribution ^{V;t) is 



k=0 



(v - V)^+^p:^^fiv; V- t)dv - / {v- Vf+-'p+^,f{v- V- t)dv 



\fe+2 + 



(39) 



(40) 



3. THERMODYNAMIC LIMIT FOR THE PISTON 



The thermodynamic hmit for the piston is defined hy A ^ oo, M oo, oo with 7 = 2mA/{M + m) and 

= N^m/M fixed. In this hmit, a = and the equations for collisions (|l]), (2) are simply: 

v'^2V-v V'^V (41) 

Similarly, the boundary conditions for the one-particle correlation functions write: 

p-{X, V- X, V- 1) = p- {X, {2V -v);X, V; t) (42) 

p+{X,V-X,w-t) = p+[X,V-X,(2V ~w)-t) (43) 

Assuming that we can permute a — > with the sum over k in Eq. (^o|), which we expect to be verified for initial 
conditions such that < F >t is not too large for all i, the equation for the evolution of ^{V-.t) reduces to: 

dMV;t) = -j^(^JJ{v~V)'p:^^f{v;V;t)dv- £ (v - V)^ pt^,f{v;V;t)dv^ (44) 

With the initial condition we have considered, it is natural to expect that it is possible to write: 

ptrfi^;V;t)^a^{v;V-,mV;t) (45) 

where {v; V; t) are non-negative continuous functions of V ^ while ^{V\ t) is a distribution in V for all t. As we shall 
see, the assumption (Esf) will be justified a fortiori after the discussion. 



Under the assumption (p5|), the equation for the evolution of (^{V]t) takes the form: 

d 



where 

F(y;t)=7 
is continuous in V for all fixed t. 



dMV-,t) = ~ — [<^{V-t)F{V-t)] 



(v-Vfa-{v;V-t)dv~ / (v - Vfa+ {v;V ■,t)dv 



(46) 



(47) 
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Property 1: 



Let V{t) — TtVo be the solution of 

d 



ViO) = Vo 



(48) 



then the solution of (pq) is 



^V;t)^^T_tV;t^O) ^{T.tV) (49) 

Proof: 

For t = 0, ( ^9| ) is an identity. For t > 0, ( ^ ) means that for any test- function g(F): 

<$t;.g>= y" t).g(T/)dl^ = ^ 0)5(nl/)dF (50) 
From Eqs. (|| - |5^) follows that: 



/ $(y; 0) (^5) (1^' = rtV)FinV; t) dV 
mV';t){^g) {V')F{V';t) dV 
-I^mV';t)FiV';t)] giV')dV' 



(51) 



I.e. 



±^V,t)^~-^MV;t)F{V;t)] (52) 

Therefore the distribution ^{V; t) defined by ( ^ ) is solution of Eq. (^). Let us remark that the continuity of F{V\ t) 
in V is essential in order that ^{V;t)F{y;t) is a well-defined distribution. 
From Property 1, we obtain immediately the following corollary. 

Corollary 1: 



If 

^{V]t = Q) ^ 5{V -Vo) (53) 

then 

^{V]t) ^ 6{V -V{t)) (54) 
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where V{t) is the solution of (|8|). Moreover, from Eq. ( ^5| ) we have 

Pturfiv; V; t) = p± ,^(«; t)5{V ~ V{t)) (55) 

where 

ptrM,t)^a^{v-V{t)-t) (56) 

and 

j PtrM^t)dv = Ptrfit) (57) 

It follows from the corollary that for the initial condition: 

'^{X,V]t^Q):^6{X^Xf))5{V -Vq) (58) 
the evolution of the piston is deterministic in the sense that 

-^{X^V-^t) ^ 5{X - X{t))5{V -V{t)) (59) 
where X{t) is obtained by integration of V{t) with X{t = Q) = Xq. 



As for ( [45D it is natural to assume that the distributions pi^p and pp^i are absolutely continuous with respect to 
'^{X,V,t), i.e. 

Pi.p{x,v:X,V-t) =b-{x,v;X,V;t)-^{X,V-t) 

pp^i{X,V-x,v-t) =b+(X,V-x,v-t)^{X,V-t) ^""^ 

where are continuous functions in {x,X,V). We obtain from (pO|), ( pO[ ) and ( pT| ) 

pi,p(x, f ; X, t/; i) = (x, w; i)*(^, V; t) 

(61) 

PKi (X, F; i;; t) =p+{x,v- t)^{X, V; t) 

where 

p^{x,v-t) = b^{x,v-X{t),V{t)-t) (62) 
Moreover, from (p^), (p3[), (pi]), we obtain for the joint distribution for the piston and particle velocity 

PtrM^y-^i)-ptrfM5{V-V{t)) (63) 
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where 

pfurfi^-^t) = P^iX{t),v;t) ^ p^^^f{t)f^{v;t) (64) 
and f^{v;t) is the probabiUty distribution for particles at the left and right surface of the piston. 

In conclusion, in the thermodynamic limit for the piston (A/ = oo, taken as described above), the motion of the 
piston is purely deterministic and the correlation functions for one particle and the piston have the factorization 
property (^l]). 
Introducing 

/•oo fV 

F2[V;ptrf{-)]^ {v-Vfp:^,f{v-t)dv~ {v^Vfpt^^f{v;t)dv (65) 

JV J — oo 

we have thus the following equations for the evolution: 

'^=im:ptrf{-)] (66) 

(dt + vdx)p~ {x, v; t) = 5{x + L)vp~ {—L, v; t) 

+ S{x - X{t))[V{t) - v\\d{y{t) ~ v)p-{X{t),2V{t) -v;t)+ e{v - V{t))p-{X{t), v; t)] (67) 

{dt + wdy)p+{y,w;t) = ~S{y - L)wp+{L,w;t) 

- S{y - X{t))[V{t) - wMw - V{t))p+{X{t), 2V{t) -w;t) + e{V{t) - w)p+{X{t),w; t)] (68) 
It follows from (|3) (||) that: 

p~{x,v;t) ^ p^{xo,vo;t^O) (69) 
p+{y,w;t)^ p+{yo,wn;t = 0) (70) 

where {xq, vq), resp. {yo, wq), is the initial conditions necessary to reach the point (x, v) (resp. {y, w)) at time t, under 
the free evolution of the particle with elastic collisions at the boundary x = —L and x = X{t) (resp y = X{t) and 

Finally, from ( |63| ) and (|6l|), the assumptions (|4|) and (|6^) will be satisfied, so that the whole reasoning is consistent. 
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Let us summarize our results in the following statement: 
Property 2: 

In the thermodynamic limit for the piston, and for the initial condition 

V;t^O)^ S{X - Xo)SiV - Vo) 

the evolution of the piston and the gases is described by the autonomous system of equations (|66| - 16^) . The evolution 
of the piston is deterministic, i.e. 

^iX,V;t)^5iX-X{t))S{V~Vit)) (71) 
and the joint probability distributions for one particle and the piston have the factorization property 

p^{x,v;X,V;t) = p^{x,v;t)^iX,V;t). (72) 

Thermodynamic quantities on the surface of the piston 
We introduce F^{V), functionals of p^^rfi'"'^'^) defined by 



V 



F^iV)^ I {v-VYp:^,f{v-t)dv (73) 

(74) 

F^{V)^ j {v-Vfp+^^f{v-t)dv (75) 



In the thermodynamic limit, 7 simply writes 7 — 2mA/ M and the equation of motion ( pGI) is of the form 

^ A [2mF-{V) - 2mF+{V)] (76) 

where 2mF2 {V) is the force per unit area exerted by the left particles on the piston and similarly —2mF2{V) is the 
force per unit area exerted by the right particles. Since -^2 (^) ^'^d ~F2{V) are monotonically decreasing functions 
of V (for any distribution to define the pressure p, the friction coefficients X{V) as well as the 

density p and the temperature T associated with the particles which are going to hit the piston by 



P"=2/ P7urfi^-^t)dv (77) 











2 / Psurfi^-^t)dv (78) 



00 



13 



2mF^{V = 0) ^p^ksT ■ 



(79) 



2mF^{V)^p^± ^\^{V)V (80) 
where \^{V) is positive for all V . The equation ( [76| ) for the evolution of the piston is thus 

|y = ^ (p - + ) - [A-(F) + A+(F)]y (81) 
and involves the difference of pressure and a "friction force" . 



We notice that for distributions p^^rfi'^'^ which are symmetric in v, (|7^), (78) and (79) are the standard definitions 



surf 

of density, pressure and temperature. However this symmetry property is not satisfied in our case. 
Moreover in the thermodynamic limit for the piston we also have from (|67|), ( |68| ) 



d f<E+ > 



2mF+{V)V (83) 



dt\ A 

where < > /A is the (kinetic) energy per unit area of the particles in the left and in the right compartment. From 
the first principle of thermodynamics we thus have the following result. 

Property 3: 



In the thermodynamic limit the evolution of the system is "adiabatic" , i.e. there is only work and no heat involved. 
Remarks 



1. It should be stressed that the deterministic evolution (|71|), the factorization property (|72|), and the adiabatic 
property ( ^2[ - |8^ ) are only valid in the thermodynamic limit. For systems with finite mass M, these properties 
will only be approximately verified for a time interval of the order of M pOf . 



2. The system of equations ( |66| - |6q ) which we have obtained is identical to the equations derived in [ [15| using a 
different point of view. Indeed, as mentionned in the introduction, the authors have considered in p5| another 
thermodynamic limit in wich L is not fixed but tends to infinity, together with ^ L^, ^ Li^. Position is 
then scaled with y = L^^ x and time is scaled with t = L^^ t. However they also introduce an assumption on 
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the distribution functions which is equivalent to a scaling on the number of particles with N = L^^ ^ L^. 
Therefore it is not surprising that their equations in terms of the scaled variables is identical to the equations 
we have obtained above. 

4. EVOLUTION OF THE PISTON IN THE THERMODYNAMIC LIMIT (M = oo) AND FOR THE 

INFINITE CYLINDER (L = oo) 

In this section, we consider the limiting situation where the length of the cylinder is infinite and the thermodynamic 
limit for the piston has been taken (Sec. 3) 

Assuming that the initial distributions /^(v) are zero if |w| ^ [vmin,Vmax], where Vmin and Vmax depend on the 
initial temperatures and pressures, there will be no recoUision of the particles on the piston. In this case P7tir/(^'^)' 
resp will be independent of t for all v > infj V{t), resp. w < supj V{t), and thus F2{V) is independent of 

t. If the initial distributions are Maxwellian with — sufficiently small, the velocity of piston will be small and 
we also expect that F2{V) is approximately independent of t. In those cases p^{v] t) — p^{v) for w ^ 0, T^{t) = 
and p^{t) — ■ 

The evolution of the piston is thus simply given by the ordinary differential equation 

iV = (A) 2mF2iV) 

(84) 

V{0) =0 

where F2{V) is a strictly decreasing function ofV : 

2mF2iV) = {p- -p+)-^\{V)V X{V) = X+{V) + X-{V)>0 (85) 

If the initial conditions are such that = p+ then from (p4|), (|85|), we have V{t) = for all t, i.e. the piston remains 
at rest and the gases on both sides remain in equilibrium with their respective temperatures and T+i the piston 
is therefore an adiabatic wall in the thermodynamical sense. 

If the initial conditions are such that ^ p+ (i.e. p^kgT^ ^ p^ksT^), then from (p4|), (|85|), the piston will 
evolve monotonically (i.e. either X{t) increases, or decreases, for all t) to the stationary state with constant velocity 
V solution of the equation F2{V) = 0, i.e. 

p-{v){v -Vfdv ~ p+{v){v -Vfdv (86) 

V J-oo 

Moreover, the approach to the stationary state is exponentially fast with time constant tq 

To = i where X = X{V = 0) (87) 

A 
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Assuming that for k = 1 and 3, the expressions 

/o 
v''p'^{v)dv 
'OQ 

are approximately the same functions of and r+ as for the Maxwelhan distributions, we have 



ro-^ = A = - y ^ (p-VT- +pWT+) (89) 



and V is given by the solution (closest to 0) of 



kB {p~T~ - p+T+)-v d'^^ ip-Vr^ + p+Vf^) 



In conclusion, the stationary velocity ^ is a function of p'^ / p^ , T+, and m, but does not depend on the values 
M/A, /M . Furthermore the time necessary to reach the stationary velocity is characterized by tq, which from 
(p7|), ( |89| ) tends to zero when /M is very large (with fixed T^, r+). In Sec. 6, we will check the above conclusions 
(in fact the assumptions under which they have been derived) by means of numerical simulations. 

To conclude this section, we remark that in the thermodynamic limit the stationary state is stable with respect to 
small perturbations of V . We should also remark that for small V and Maxwellian distributions for the velocities, 
the equation (^6|) is similar to the equation obtained from fluid dynamics under the assumption that the evolution is 
adiabatic 0. 

5. EVOLUTION OF THE PISTON IN THE THERMODYNAMIC LIMIT {M = oo) AND FOR A FINITE 

CYLINDER (L < oo) 

If the cylinder has finite length, we have to consider the full set of equations ( p5[ - p8| ). In this section, we shall give 
only a qualitative discussion and then in Sec. 6, we shall compare our conclusions with the numerical simulations. To 
simplify the notation, as in the simulation we now place the left boundary at a; = and the right boundary at x = L, 
i.e. the cylinder has now a finite length L. 

Introducing the average temperatures of the gases in the left and right compartments by 

the conservation of energy implies in the thermodynamic limit for the piston (recall that < V'^ >=< V >^ in this 
limit) 
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where N ^ N+ + N' . From Eqs. (xm, E2, 83), we then have 



jV^I^^) 2m [F,-{V) - F+{V)] (93) 

kBjT- = -4m (^A ) F,-{V)V (94) 

fcB|r+-4m(^A^ F+(l^)l/ (95) 

where 

2mi^2'^(F) = p^ksT^ ± (^^^ A±(V^)T/ (96) 

We notice that ( |9^ ) is satisfied fo aU Let us first remark that in (|93|), appears iTi{A/M)pf^^^ which at < = is 
R-/Xo and R+/{L-Xo) with 

= m — (97) 

Therefore, as discussed in Sec. 4, if is sufficiently large then after a time tq = 0{{R^)~^) 1 the piston will 
reach the velocity V = V{p'^ / ,T^) given by Eq. (p6|), i.e. X(t) = Xq + Vt. This evolution will continue over a 
time interval of the order 

. . 2Xo 2{L - Xq) . 

n = mm ^ -, = ^ (98) 

3bT /3bT± _|_ y 



which is the time for the sound wave to propagate from the piston to the closest boundary of the cylinder. 

For t > Ti, the density distributions P^urf^'"'^^)^ w > depend explicitly on time. To discuss the evolution we now 

introduce the following 

Average assumption: 

The thermodynamic quantities associated with the surface of the piston (|77| - ^o|) , defined by pf^rfi"^'^) with u ^ 0, 
are approximately equal to the average of the corresponding quantities in the left/right compartcment, i.e 

f±(t)~r±(i) (100) 
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Let us remark that this average assumption is usually introduced in the experimental measurements of the ratio of 
the specific heats of gases, 7 = Cp/ci, [4,18]. 
With this assumption, Eqs. (pslM) yield 



d N~ T~ iV+ T+ 



(101) 



d V . f M \ \-{V) 

dt X \N- kg 



(102) 



L-X 



dt L-X \N+ J ks 

(once again one can see that (p2) is satisfied for all t). Therefore introducing 



(103) 



a- = VT-X 



and 



a+ = Vt+{L - X) 



(104) 



eqs. ( |101[jl03| ) implies 



dt 



a 



X-{V) 



AT 

A I p-keVT- 



(105) 



dt 



X+iV) 



/T+ 



(106) 



Let us note that the entropy of the "one-dimensional" perfect gases is 



^ = ) fcB ln[VT-X] + fiN-/A) 



(107) 



Therefore from (105), (106), 



S+ fN+ 



A V A 



kB in[Vr+(i - X)] + /(iv+M) 



d_ fs^ 
dt[ A 



M\ A±(F) ^2 



A / r± 



and we have thus obtained the following 



Property J^: 



(108) 



(109) 



Under the average assumption ( |99[]lO0| ), the evolution of the piston is adiabatic, i.e. no heat transfer is involved 
and the entropy of both gases defined by (|l07|),(|l08|) are strictly increasing in time. 
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Moreover, assuming that 

/ vPsurfiV'i)dv and / v''pf^^f{v;t)dv 

Jo J-oo 

are approximately the same functions of and r+ as for the MaxweUian distributions, we have (89) 



I.e. 



N- 



TT X2 



A+(y = 0) 



M 



TT {L-xy 



In conclusion, from (105), (106) and (111) we have 



^(a--c7+)=0(F^) 



N-\ [a-f 



[a+f 



M J \ M J {L- Xf 

and thus if the evolution is such that the velocity is small for all i, then 

(7^ — (7^ = cte 



\{V)V 



(110) 



(111) 



(112) 



(113) 



(114) 
(115) 



(116) 



Equilibrium point (M = oo): 



From Eqs. (^,101-103) we find that the equilibrium point is given by 



N- 



N I f 



T7 = To 



L 



i^;^/ = ^o(i- 



L 



i.e. in particular p^. = pjr . Furthermore from (114) taking a — w cte yields 



IT^X - Vf^{L -X)=C= ^T-{0)Xo - ^T+{0){L - Xo) 



(117) 
(118) 



(119) 



i.e. with (117), (118) 



N 



m)(.-.v,3.;|c 



(120) 



Solving (120), (117), ( |llq ) with the constants C, Tq given by the initial conditions (|92[), gives the equilibrium 
state {X f jTjT ,T^) which is a state of mechanical equilibrium pj = p'j^ but not thermal equilibrium Ty ^ . 
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It is important to realize that thermostatics yields only the condition of equality of the pressures, i.e. mechanical 



equilibrium. In our dynamical approach, we have obtained a new equation (114), consistent with the second principle 
of thermodynamics [property 4 ) which then enables us to determine the equilibrium point. 



Linearization around the equilibrium point: 



Let {Xf, Tf) be the equilibrium point given by (117-120) and x = X — Xf. Linearizing the equations (101), (105) 



(|106D, together with (111) for the friction coefhcient, yields 



with 



fT ~ cte 



(7+ = cte 



TJ{L-Xs) 



(121) 



X = —UJqX — \x 



N 



ksTo 



MJ Xf{L-Xf 



(122) 



(123) 



A 



ld,m N ksTo 
TT V M L 



IN- 1 

Iw" X] 



lN+ 



1 



(124) 



M {L~Xf) 

In conclusion the equilibrium point (for the thermodynamic limit) is stable and the approach to equilibrium is a 
damped, harmonic, oscillation 

A 



X{t) = C+e^^ ^ + C^e 



(125) 



Moreover, in the linear approximation the evolution of the gases are at constant entropy. 



Remarks 

1 . Although our derivation of the evolution equation was strictly Hamiltonian (Liouville equation) , and the concept 



of entropy was never introduced, it is interesting to note that the frequency cjq, (122), which we have obtained 
with our assumption coincides with the frequency obtained in thermodynamics assuming "adiabatic oscillations" 
[4,18]. 

2. For the case — = N which will be considered in the simulations (Sec. 6) 

^=^ = ^ (126) 
MM ^ ' 
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To = -{T-{0)+T+iO)) 



(127) 



8R 



L \ L 



(128) 



This implies that the motion is weakly damped if 



with "period" 



27r 




L 



(129) 



Wo ^1 - 



(130) 



We have thus obtained the conclusion that if i? < 2.3 the motion is weakly damped, and if i? > 4.7 the motion 
is strongly damped. Again we should mention that the main result of the recent experimental measurements of 7 
is that one should distinguish between two regimes, corresponding to weak and strong damping with very different 
properties [18]. In particular, experimental results show that the frequency of oscillations for weak damping is very 
close to the values obtained assuming adiabatic oscillations, eq.(122). 

6. NUMERICAL SIMULATIONS 



To illustrate the results of Sec. 4 and to check the conclusions of Sec. 5 which were based on the average assumptions 



(|99|), (100), we have performed a large number of numerical simulations. In all simulations, we have considered a 
one-dimensional system of fixed length L, with the left boundary at a; = 0, and we have taken: 



fcs = 1, m = 1, L = 60- 10^, Xa = 10 . 10*, Vq = 



(131) 



N+ = N, 



R 



N 



p-(0) = 5p+(0) 



(132) 



A very large number of simulations has been conducted with the number N of particles in the left /right compartments 
ranging from 10'^ to 500 • 10'^, the mass M of the piston from 1 to 10^, and the parameter R from 0.1 to 10^. 

At the initial time t = 0, the left and the right particles are taken with Poisson distribution for the position and 
Maxwellian distribution for the velocities, characterized by: 



T" = 1, T+ = 10, i.e. To 5.5 



(133) 
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and thus 



p+ = 2p- 



(134) 



Since under the scahng V = aL, Xq = aXo, time is scaled with t' — at, we have kept L fixed. Similarly under the 
scaling T'* — aT^, time is scaled with t' = ^^^d thus we have also kept T~ and r+ fixed. 

From the discussion of Section 4, it is expected that for R sufficiently large, the actual velocity of the piston before 
the first recoUisions of the particles will coincide with the stationary velocity V, given by eq.(pq), for the infinite 



cylinder and the initial conditions (133-134). For the parameters used in the simulations eq.(86) yields 



V = -0.3433 



(135) 



The velocity c of the sound wave in the one-dimensional perfect gas is given by c = ■\/3fcBT/m. Therefore for our 
simulations, the time needed for the sound wave to make the first collision on the piston is: 



= 0.96 . 10*^ 



c- ~V 

The mechanical equilibrium defined by eqs. (117 1201) is: 



t+ = ^ 1.95 . 10^ 

c+ + V 



(136) 



Xf = 8.42 . 10^ 



T- ^ 1.54 



T+ = 9.46 



p- = 1.83(4) N . 10" 



(137) 
(138) 
(139) 



For the approach to equilibrium, we have from eqs. (122), (128): 



Wo 



\fR . 0.2756 . 10" 



R. 0.331 . 10" 



(140) 



and the "period" of the damped oscillations is: 



T = 2.28 . 10^ [R{1 - 0.36i?)]-i/2 



(141) 



which yields: 



R„ 



(142) 



Let us remark that for thermal equilibrium of a heat-conducting piston, one would have: 



Xth = 30 . 10*, T-^ = T+ = 5.5, p- =p+ = 1.833 N . 10"^ 



(143) 
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In figure 1, we have taken i? = 4 and investigated the thermodynamic Umit by considering increasing values of M 
from 10 000 to 100 000. From figure 1, we conclude that for t < 30 . 10^, the thermodynamic limit for the piston is 
reached when M > 50 000. Other simulations with fixed R ranging from 0.1 to 300 confirm this result: for t < 30.10^, 
the evolution corresponds to the thermodynamic limit if M > 50 000. 

In figures 2 and 3, we have considered the time evolution for the piston for different values of R, and values of M 
such that the evolution corresponds to the thermodynamic limit. The simulations {Fig.2) show that the evolution is 
very weakly damped for R < 1 and strongly damped for i? > 4 in agreement with the conclusion of Sec. 5. Moreover 
the piston evolves toward the equilibrium position Xf = 8.35 • 10"*, which is in excellent agreement with the predicted 
value Xf = 8.42 • 10* obtained from the result of section 5, i.e. eq.(140). 

The weakly damped oscillations have a period which is in very good agreement with the expected values of Sec. 5 
(see Tablel), but the damping coefficient computed over more than 1000 oscillations gives a value several order of 
magnitude smaller than A given by eq.( |l4C| ). For example for R — 0.2, M — 100 000, we observe a damping equal 
to 3.10^'''A. The oscillations occur around the equilibrium position Xf, and their period remains constant in time (at 
least up to i = 7000 • 10^ which corresponds to 1500 oscillations and is the time the simulation has been running). As 
predicted, the frequency of the oscillations increases with R, but is independent of M for sufficiently large M. The 
fluctuations on the amplitude of these oscillations is of order 1%. 

On the other hand for i? > 10 and M > 6 000, i.e. for strongly damped evolution, it is seen on figure 3 that the 
evolution is independent of R and M ii t < 20. 10^. We also notice (Fig. 3) that for R > 10, the piston acquires 
almost immediately a velocity: 

Vi = -0.34 (144) 



which is in perfect agreement with V computed above with eq. (135). Therefore for strong damping the friction 
coefficient coincide with the computed value during the initial evolution, i.e. as long as the influence of the boundaries 
(x — 0, X = L) did not appear. The piston then arrives at its first minimum position at the time ti — 0.9 . 10^, 



which corresponds to the time t~ for the sound wave to return on the piston, eq. (136), and reaches the position 
of mechanical equilibrium after two oscillations, i.e. at a time t which is about Ati (see also Fig. 4). It is seen on 
figure 1 and 2 that for R between 4 and 10 new oscillations appear after a time of the order 20 ti with amplitude 
SX/L sa 4 . 10^"^ and period approximately equal to 2ti independent of R and M (for AI large). The amplitude of 
these new oscillations is modulated with a period about 20 ti and they tend to disappear with time. We are led to 
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conjecture that these oscillations are associated with sound waves propagating in the gases. As we have checked, these 
oscillations reflect the fact that the velocity distributions of the gases are not Maxwellian. They will be responsible 
for the approach to Maxwellian distributions. Similar results and conclusions have been obtained in [19]. For larger 
values of R (e.g. i? = 50 or 100 as shown in figure 3) the number of particles necessary to observe the adiabatic 
evolution over a time interval 40 ti is several millions. These simulations would need a very large computer time and 
have not been conducted. 

In figure 4, we show the time evolution for the temperature of the gas and for the pressure defined by 

_ NkBT-{t) + NksT+jt) 

P = X{t) P = {L-X)(t) (^4^) 

where = 2E^ /Nks, being the energies of the two compartments. It is seen that these evolutions follow the 
evolution of the piston. Moreover, we notice that the temperature increases (resp. decreases) under compression 
(resp. expansion) of the gas, which reflects the fact that the evolution is adiabatic. We recall that the area A does 
not play any role. From figure 4, we obtain the numerical values for temperatures and pressure when mechanical 
equilibrium is reached, 

T7 = 1.52, T+ = 9.48, p- = p+ = 5.5 (146) 



which are again in perfect agreement with the expected values (137 139) — 1.54, = 9.46, pf — 5.5 since 
N = 300 000. 



Finally from (107-108), we have for the change in entropy between the initial and final states: 



AS*- = NkB 23 . 10-3 (147) 
AS+ = NkB 6.7 . 10-^ (148) 

We thus observe numerically an increase in entropy for both gases which once more confirms the adiabatic property of 
the piston. The result observed in figure 3 that for i? > 10 the evolution is independent of R leads us to the following 

Conjecture 

In the thermodynamic limit, the evolution for i? — > oo is given by the solution of: 

{dt+vd^,) p-{x,v;t) = 6{x) V p-{0,v;t)+6{x-X{t)) [V{t)-v] p- {X {t),v;t) 
(dt+vdy) p+iy,v;t) = ~Siy - L) V p+iL,v;t) ~ S{y - X{t)) [V{t)~v] p+iX{t),v;t) 
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Jv(t) 




This conjecture follows from (|65|-|68|), with the left boundary placed at the origin (a; = 0), and the boundary condition 



at a; = X{t) 



p^{X{t),v;t)^p^{X{t),2V{t)-v;t) 



for V < V{t) 



which is expected to hold in the thermodynamic limit IT 



7.CONCLUSIONS 



As mentionned in the introduction, we have shown that in the thermodynamic limit for the piston [A oo with 
N^/A, M/A, and L fixed) the evolution of the piston is deterministic, i.e. the distribution function is 'i>{X,V;t) — 
6{X — Xt)d{V — Vt) where X{t) is the solution of a system of autonomous equations coupled to the one-particle 
distributions for the gases. Moreover we have shown that in this limit the two-point correlation function for one 
left (right) particle with the piston factorizes as the product of the individual distribution functions. Furthermore, 
the evolution is strictly adiabatic (for M — oo) in the sense that the changes in energy are entirely due to work, no 
heat transfer is involved and the entropy of both gases are strictly increasing. If the length of the cylinder is infinite 
(and M — oo), the piston evolves towards a stationary equilibrium state with velocity V approximately proportional 
to (p+ — p^). Therefore, in the thermodynamic limit (and L = oo) the stationary state is a state of mechanical 
equilibrium, iff j5+ = . 

If the length of the cylinder is finite, we have introduced an assumption to express the density of particles and the 
temperature at the surfaces of the piston by their average values in the left and right compartments. With this 
assumption, we were able to analyze the evolution and the final state. The numerical simulations, as well as the 
analytical expressions, have shown that for finite L the motion is characterized by damped oscillations. It is weakly 
damped if N/M is small and strongly damped if N/M is large. Moreover the numerical values obtained in the 
simulations confirm all the conclusions derived from the average assumption, except for the damping oscillations 
which appears a lot smaller than predicted. We are thus led to the conjecture that there are two different mechanisms 
responsible for friction. One mechanism correctly described by the discussion of Section 4 is associated with the 
motion of the piston in the absence of recoUision {L = oo) and is responsible for the stationary state in the infinite 
cylinder. Another mechanism which is responsible for the damping of oscillations appears to be associated with sound 
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waves, i.e. inhomogeneities in the gases, traveling back and forth between the piston and boundaries. Finally we have 
observed that over a time interval of order M the evolution is independent of N/M, as soon as N/AI is large enough. 
It remains however to give a proof of the result without introducing the average assumption. 

In a forthcoming paper, we shall present numerical simulations which show the key difference between the infinite- 
mass and the finite-mass problems pc[ |. For the finite-mass case, we have the following picture. In a first stage, 
characterised by a time scale of order ti = L^/mJEq, the motion of the piston is adiabatic (no heat transfer) and 
corresponds to the motion that we have described in the thermodynamic limit: it is deterministic, with either weak or 
strong damping, temperature increases under compression; the evolution is independent of M (for M large enough) 
and it proceeds until mechanical equilibrium j<+ = is reached. In the second stage, which is now characterised by a 
time scale of order T2 — Mri/m, the evolution has exactly the opposite properties: it is stochastic, which implies heat 
transfer through the piston, and the system evolves with constant pressure i.e. p^{t) — Pf + 0{1/M), to a state of 
thermal equilibrium {T^ = T+); in this second stage the evolution depends strongly on M , more precisely we observe 
the scaling relation XM{t) — X{t/M) where X{t) is independent of and M. 
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FIGURE CAPTIONS 



Figure 1: Thermodynamic limit. Position of the piston in function of time fixed ior R = N/M = 4 and 

a) M = 10 000, b) M = 25 000, c) M = 50 000, d) M = 100 000. Predicted mechanical equiUbrium for adiabatic 

piston Xf = 8.42 10^. Thermal equilibrium for conducting piston Xth = 30 10^. 

Figure 2: Prom weak to strong damping. Position of the piston in function of time, a) i? = 0.1 and M = 100 000 , 
h) R = 0.2 and M = 100 000 , c) i? = 4 and M = 100 000 , d) i? = 10 and M = 30 000 . 

Figure 3: Strong damping. Position of the piston in function of time, a) i? = 50 and M = 6000 , h)R = 50 and 
M = 10 000, c)R = 100 and M = 4000 , d)R = 100 and M = 6000. 

Figure 4: Adiabatic evolution: approach to mechanical equilibrium for i? = 10 and M = 30 000. (a) Position, (b) 
Pressure, (c) Temperature. The predicted value for the mechanical equilibrium are Xf = 8.42 10*, pj = P~f = 5.5, 
= 1.54, T+ = 9.46. 

Table 1: From weak to strong damping. Simulations with M sufficiently large to describe the thermodynamic limit. 
Value of the first minimum of the position of the piston and time at which it occurs. Comparison between the observed 
and predicted period of oscillations. 
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